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Abstract 

We discuss how to apply the dimer method to Ising models on group lattices 
having non trivial topological genus g. We find that the use of group extension 
and the existence of both external and internal group isomorphisms greatly 
reduces the number of distinct Pfaffians and leads to explicit topological for- 
mulas for their sign and weight in the expansion of the partition function. 
The complete solution for the Ising model on the Klein lattice group L(2, 7) 

with g = 3 is given. 
PACS N. 05.50 
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I. INTRODUCTION 



Among the known approaches to the exact evaluation of the 2D Ising partition func- 
tion which have followed the celebrated Onsager solution jl]], the dimer method fully 
exploits the combinatorial and group theoretical properties of the lattices by relating the 
partition function Z to the generating function Z^ of close-packed dimer configurations. 
Though the method is in principle independent from the dimensionality of the lattice, the 
corresponding analysis in three dimensions has never been developed, due to the difficul- 
ties in extending Kasteleyn's Theorem on lattice orientation. In this paper, we deal with 
the issue of generalizing the dimer approach to the case of non-abelian 2D lattices of high 
topological genus which should in fact be equivalent to higher dimensional lattices. Indeed 
the 3D cubic lattice can be considered as a handlebody 2D lattice of genus g = N/4 where 
N is the number of sites of the lattice. This hints that a non vanishing ratio g/N may be 
related to an effective dimension D > 2 of the lattice. The great difficulty of the problem 
suggests that a possible concrete way to analyze such lattices is to consider at first graphs 
possessing the largest possible symmetry group, and the best candidates appear to be the 
group graphs already well-known in the mathematical literature In this framework 

we may consider the 2D planar lattice with periodic boundary conditions as an example of 
abelian lattice group (translational symmetry group) with genus g = 1, whereas an example 
of finite non-abelian lattice with g = is given by the fullerene-like lattice studied in ref. 

The paper is organized as follows. In Sec. II we outline the basic ideas concerning ap- 
plication of the dimer method to group lattices. In Sec. Ill we define an extended lattice 
group H and relate it to the homology and cohomology groups (mod2) Hi, H 1 of the original 
lattice A, define the Grassmann algebra over the decorated lattice A* and the Pfaffians as 
function on orbits of H 1 . In Sec. IV we apply the results of Sec. Ill to discrete groups of 
particular interest such as the Klein group L(2, 7) (of order A" = 168 and genus g — 3) 
and discuss the role of external automorphisms. In Sec.V we analyze in detail the orbits of 
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Hi,!! 1 under the action of H, define an invariant duality map </? : Hi — > H 1 and auxiliary 
functions of relevant computational interest. Finally, in Sec. VI we construct explicitly the 
irreps of H, apply them to the harmonic analysis on A* and derive the expansion for Z. 
Few preliminary numerical results are also given. 

II. GROUP LATTICES AS ISING LATTICES 

We recall here key points of the dimer procedure which are specific to group lattices (see 
also refs. g§) : 

(1) . A discrete group G is defined by a presentation given in terms of a set of p generators 
Ak, k — 1, p and t relators p i — 1, . . . , t. The P, are words in the generators equivalent 
to the identity. Let T be the free group on the k = 1, . . . ,p. Let Af be the minimal 
normal subgroup of T containing all relators. Then by definition G = T jH ■ 

(2) . The Cayley lattice A for a group is defined by giving a map L : G — > R 3 where 
L(g) G R 3 is the point of R 3 corresponding to the group element g G G. A bond of color 
k is then a line joining L(g) to L(Akg). The genus g of A is that of the surface of minimal 
genus S C R 3 on which A can be drawn. The Ising Hamiltonian is then defined as 

^ p 

E = ~ H J k Oh<?A k h , (1) 

h&Gk=l 

where the {<Jh = ±1} are the spin variables and {Jk} the exchange interactions between 
connected spins. 

(3) . Each relator Pj is then represented on A as a closed circuit ((Pi) made of oriented 
colored bonds. If C(P) encloses a simply connected region (tile) on S then Pj is called local 
relator. If S has genus g = then all relators are local. Of particular interest are models 
where g is large. 

(4) . The group lattice A is interesting on its own but cannot be related directly to the 
partition function of an Ising model and to do this we must consider Q a decorated lattice 
A*. This amounts to replace each site of coordination q (q > 2) of the original lattice by 
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a sublattice containing 3(q — 2) points and Aq — 9 decorating bonds. The Ising partition 
function is then related to the dimer covering generating function on the decorated lattice. 

(5) . In order to compute the generating function we orient A* according to the Kasteleyn 
prescription by assigning arrows to each bond inherited from A in such a way that for any 
closed circuit £ on A # , the number of bonds of £ oriented clockwise is of opposite parity to 
the number of sites enclosed by I. For the decorating bonds see || or Sec.V. 

(6) . The Kasteleyn rules define completely the orientation for lattices of genus g = 0, 
whereas for lattices of higher genus we have to deal with further sign fixing for loops which 
are not homologically trivial, i.e. not the boundary of a union of tiles. The assignment of 
arrows to A* (or A) is not invariant under the action of G but rather under an extension H 
of G closely related to the homology H% and cohomology H 1 groups of A. 

(7) . The dimer covering generating function of the lattice can be expressed as a weighted 
sum of Pfaffians Pf{4>) , where <fi G H 1 and with sign given explicitly by the if -invariant 
function 6((), ( £ Hi defined in Sec.V. Harmonic analysis on H allows to factorize Pfaffians 
into determinants of lesser order and external automorphisms induce identifications between 
Pfaffians. 

III. THE EXTENDED LATTICE GROUP 

In this section we first discuss the group extension H of G. Next we show that H 
partitiones the homology H\ and the cohomology H 1 groups of A into non-intersecting 
orbits characterized in terms of sign functionals. Their role in the expansion of the dimer 
generating function is then analyzed. 

A. The groups H, H\ and H 1 

The extended lattice group H can be obtained from G by replacing the relators P{ with 
new relators containing the following elements: 
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(i) If Pi,Pj G G are local, then P i P j 1 , Pf, P i A k P i l A k 1 are relators in H; 
(ii) For generic relators P i: Pj G G then Pf, Pj, PiPjP^Pj 1 are relators in H. 

We can write then Pj = Q (Q 2 = 1) for all local relators Pj, with A k Q = QA k , i.e. Q 
is a central element that we call central signature. Putting Z { = P { Q the non local relators 
will be written as Zf, ZiZjZ~ x Zj 1 where in general A k Zi ^ ZiA k . The {Zi} generate an 
abelian normal subgroup HZ C H . Particular examples of this extension will be discussed 
further on. 

Non trivial loops form a chain group Ci(A # , Z 2 ) . The equivalence classes of Ci(A # , Z 2 ) 
modulo boundaries form the homology group Hi(A#, Z 2 ). The class of multiplicative func- 
tional on iJi(A # ,Z 2 ) with values ±1 are then the elements of the cohomology group 
H l (A*,Z 2 ). Notice that H l (A*,Z 2 ) ~ H l {A,Z 2 ) and H 1 (A*,Z 2 ) ~ H 1 (A : Z 2 ) and there- 
fore we denote them briefly by the symbols H 1 , Hi respectively. When dealing with elements 
of HZ, H 1 we may use addition instead of the product as composition rule We write g ~ g' 
if did' £ H define the same site on A. 

H\ is isomorphic to HZ . To see it consider the closed chain t on A as defined by the 
sequence g = g, g p — h p g p _i, p = 1, . . . , n with h n+p = h p . is then defined as 

CW = 9Q 1 KK-i---h 1 g = g^hp^hp^. . .h 1 h n h n _ 1 ..h p+1 g p . (2) 

Starting from g we move n steps across bonds in A and each bond defines an element 
hi G H, where hi = A^ 1 , A k being a generator of H. If g n ~ go, £ is closed, ((£) G HZ 
does not depend on the choice of g on £ and defines a map ( : Ci(A, Z 2 ) — > iJZ. Adding 
a boundary £ to £ (i.e. adding the boundary of a union of tiles on the lattice) amounts to 
replace a sequence of bonds h^ , . . . , by an equivalent one obtained by using local relators 
only , therefore ((£ + £$) = ((£) and ( induces the map k : Hi — > HZ. Given inversely 
an element ( G HZ, expressed in terms of relators as £ = h n . . . h±, a corresponding closed 
chain G Ci(A, Z 2 ), and therefore a cycle in iJi(A, Z 2 ), is given by a sequence of gi, 
i = 0, . . . , n, with g = 1 and = Also if C = ^ n • • • h, (' = /i m +n • • • h +n G HZ 
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and g n ~ g m+n ~ lthe chain £(( + (') = £(() + £((') G C 1 (A,Z 2 ) denned by g = land 
hiQi-i — gi,i — l...m + n corresponds to the element £ + £' G iJZ. Hence k is a group 
isomorphism. 

Right multiplication on A by an element h G H translates the group lattice and if we 
replace with g^li in i we obtain a closed path £/j which is the right translation of I by h 
and a corresponding element C(-^i) = h~ l I h which defines the action of H on I G HZ. In 
this way H acts naturally on A* , H\ and if 1 and partitions H x ,Hi into non-intersecting 
orbits. The use of orbits greatly simplifies the computation of Z. 

B. Sign functionals and lattice orientation 

We assign an orientation to A* according to the Kasteleyn rules and to each site h G H a 
Grassman variable a(h), with an anticommuting wedge product a(h) Aa(h') = — a(h) Aa(h'). 
Reversing the arrows on all bonds sharing the same site i corresponds to the change a(h) — > 
— a{h). 

Let I given by (0). Each hi is of the form A p ^. where pi = ±1 determines the orientation 
of the arrow in the bond. Closing I means that we must identify a(g ) and a(g n ) as 

a(9o) = P n a(g n ) (3) 

We define then 

n 

v(C(A) = -UPi ■ ( 4 ) 
i=i 

If £ is trivial then r](() = 1. All rj(() defined in this way, hereafter called sign functionals , are 
characterized by a particular recursion relation which can be proved as follows. Following 
ref. we introduce a reference dimer configuration Co on A* consisting of all bonds of A* 
inherited from A. By definition, any other dimer configuration C when superposed on C 
generates transition cycles £ G H\. In this proof we use transition cycles only. We consider 
representative chains £,£',£" G Ci(A # ,Z 2 ) of the pair of intersecting cycles (,(' an d of the 
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cycle ( + (' as in Fig.l. Chains £, £' run over the sequence of sites Pi, i = 0,1,2,7, 8,9 and 
i =5,6, 7, 2, 3, 4 respectively, while £" runs over % — 0, 1, 3, 4 and i — 5, 6, 8, 9 but not i = 2,7. 
We denote with the bonds joining Pi, Pj and with n (£) the number of anticlockwise 
arrows on £. The key point is that bonds B\ 2 , -B13 are not contained in the intersection of 
( and C + C an d thus Bi 2 , B ri ^ C whereas -B i G Co- Similarly B 3i , B 27 , B 5G , B 89 e C . 
By orienting A* according to the Kasteleyn rules and counting arrows we see that n (£") = 
n(£) + n (£') + 1 (mod2). The general case where Bij are replaced by sequences of bonds 
corresponds to the addition of boundaries to £, £', £" and leads to the same formula. Hence 
n(£) is a function in Hi and ^(C + CO = ^(OviC) if C> C intersect and f/(C + C') = viOviC) 
otherwise , i.e. 

v(C + C) = v(CHC)(-ir {u ' ] , (5) 

where Q ((, (') = 0, 1 is the intersection number (mod2) of (, Given a sign functional 
7](Q then <f)(()i](() , <p G H 1 , is another sign functional and so is the if -invariant function 
6(C), (eHx defined in Sec.V. 

C. Pfaffians and dimer covering generating function 

Consider now the form 

/ = - J2 x hh , a(h) A a(h') , (6) 

h,h'£A# 

where Xhh' = if h and h! are not connected in A* (neighbors). The factor | is inserted in 
order to avoid double counting, the activity is given by Xb = Xhh' = —%h'h = coih(f3Jhh') if 
the bond b is oriented with the arrow from h to h', (5 is the inverse temperature and Jhh' 
the exchange interaction between spins ( in the previous notation) which depends only on 
the color k of the bond. Under the stated conditions decorating bonds have Xhh 1 = 1- For 
the sole purpose of safe handling of signs we assume initially Xb > . All arrows on A* and 
all signs in @ are then fixed by the sign functional rj. The Pfaffian Pf{rj) is defined by 
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jM/2 _ 2 M M\Pf(rj) a(l) A a(2) A ... A a(M) . (7) 

where M = 3(g — 2)iV, iV is the order of the group (or lattice size) and q is the lattice 
coordination. From ([7p one derives the well-known relation 

Pf( V ) 2 = Det(X(rj)) (8) 

where X(rj) is the M x M matrix of elements x^. In general Pf(rj) is multilinear function 
of the Xb considered as independent variables. The partition function Z is then given by 

Z = Z d 2 N X\sinh((3J a ) l l 2 , (9) 

a 

where a runs on all the oriented bonds of the undecorated group lattice and 

Z d = 2~ 9 J2 s v P f(v) 8 V = ±1 . (10) 
v 

To each dimer configuration C we associate the contribution of C to Pf(i]), a monomial 
M(C) = fib x b where b runs over all oriented bonds of C. Because of our conventions M(C) 
has sign rj(Q where ( is the superposition of C and Co- I n general ( is the sum of non 
intersecting cycles C — J2id so that 77(C) = YliT](Q). The signs s v must then be chosen in 
such a way as to set equal to 1 in ( |TD| ) the coefficient of all M(C) appearing in Once 
this is done in Sec.V the Xf, can be given any sign. 

IV. 2D ISING LATTICE AND THE £(2,7) LATTICE GROUP 

The above group extension procedure can be naturally applyed to a wide class of lattice 
groups. Here we discuss two examples: the abelian 2D Ising lattice (of genus g = 1) and the 
L(2, 7) Klein group. The latter is non-abelian and of genus g = 3, and its analysis should 
hopefully be of interest in the analysis of more general structures of dimension D > 2. 

A. 2D Ising lattice 

The Onsager solution for the 2D Ising lattice made use of a rectangular n x m lattice 
with sites labelled by i = 1, . . . , n and j — 1, . . . , m ( n even). In this lattice we identify 
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opposite sides, giving it a toroidal ( g = 1) topology and turning in a group lattice G nm . This 
last property makes it possible to apply harmonic analysis, i.e. Fourier transform methods, 
which eventually lead to the final formula. The group G nm of the lattice is defined by the 
presentation 

ST = TS, S n = l, T m = l , (11) 

for all integers n, m. However in order to satisfy the Kasteleyn rules we must use a central 
extension H nm of (|TT| ) defined by: 

ST = QTS, S n = Z 1 , T m = Z 2 , 
Q 2 = Z\ = Zl = l 
QZ 1 = Z 1 Q, QZ 2 = Z 2 Q, Z X Z 2 = Z 2 Z X . (12) 

Sites are labeled by elements h G H nm with the condition 

a (hQ) = -a(h) . (13) 

The H nm lattice is actually a covering of the G nm lattice with corresponding variables a(h) 
identified by the above condition which embodies the Kasteleyn rules. 
Besides this conditions we must set 

a(h Zi) = ei a(h) , a(h Z 2 ) = e 2 a(h) , (14) 

where €\ and e 2 take the values ±1 in all the possible combinations. These correspond to 
the 2 29 = 4 possible choices of G H l . It is tempting but misleading to label spin sites on 
the lattice by elements in G nm . In this case the orientation would no longer be invariant 
under G nm , a clear sign that the true symmetry is that of the extended group H nm . 

B. Non abelian group lattices of genus g > 1 

We examine now other types of symmetries which require non-central group extensions 
and which are of interest for models in D > 2. The natural further step is provided by 
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a vast array of discrete groups many of which are discussed in detail in the literature (see 
|5|-[8|] ). Of particular interest are the Klein groups L(2,p), where p is prime. The non- 
abelian group of genus g = (the fullerene lattice) studied in ref. corresponds to p — 
5. The next interesting example which is at the same time non abelian and has a non 
trivial topology is given by the lattice group L(2, 7) (also called T(2, 3, 7) in the context of 
hyperbolic tessellations). The L(2, 7) group , briefly called G, has order 168 and is defined 
by the presentation 

U 7 = 1 , V 2 = 1 , {UVf = 1 , (VU 3 ) 4 = 1 . (15) 

The analysis of a lattice possessing this kind of symmetry is one of the chief results of this 
paper which hopefully opens the way to the investigation of more general and interesting 
structures. The group lattice given by (|1^) has genus g = 3 (produced by the non local 
relator (VU 3 ) 4 ) and can be tessellated by 24 heptagons and 56 hexagons, for a total of 168 
bonds of type U and 84 bonds of type V (see ref. |L1|, pag. 539-549, for more details on the 
lattice). 

The group G nm given by (|IT]) is abelian and has a non trivial genus g — 1. On a lattice 
of genus g we expect that the Kasteleyn rule determines the orientation up to 2g signs with 
a total of 2 2g configurations. Each close path on the lattice A can be considered as an 
element of the homology group Hi and an assignment of all remaining signs as an element 
of the cohomology group H 1 . It is therefore sufficient to assign signs on a suitable basis of 
2g elements of Hi. In the toroidal case this was done by giving ei , e^- in fact Z\, Z 2 are 
central in the group G nm . 

The group G, as defined by ([IB]) , is non abelian and has a non trivial genus, in particular 
the Zi are no longer central elements and we deal with the discussed non central extension 
H of G defined by : 

U 7 = Q, V 2 = Q, {UVf = Q , 

W = Q(VU 3 ) 4 = Z(l) , (16) 
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and additional relators which can be expressed in term of the auxiliary elements 

Z(h) = hT 1 Wh, heH, (17) 

characterized by 

Z(h 1 )Z(h 2 ) = Z(h 2 )Z(h 1 ), (Z(h)) 2 = l, \/hi,h 2 ,heH. (18) 

The Z(h) are not all independent and can be expressed in terms of the subset of 2g = 6 
elements 

Z n = Z(U n ) = U~ n Z(l) U n . (19) 

Clearly Z n = Z ni ( mo d7) and Z = Z(l) = Z 7 . The following identities equivalent under ( |TS"D 



Z$Z\Z \ZqZ%ZqZ2 — 1 , Z\Z^Z^ZqZ2Z \Z§ — 1 , (20) 
reduces to 6 the number of independent elements Zj and can be proved by repeated appli- 



cation of relators (1T0) but not of 



A generic £ G H\ can be always written in the additive form J^i^Zi where rii = 0, 1. 
Because of J2i=o Z{ = we can always choose the rit in such a way as to have J2i n i — or 
1. All elements Z(h), h 6 if, can be obtained by repeated conjugation of the Z n by U and 
V. Conjugation by U is trivial, i.e. f/ _1 Z n U = Z n+ x. Conjugation by V is less obvious. By 
using flT6|) and not (|18|) we find that 

v~ x z Q v = z 4 , = zr 1 ^ 1 , v- 1 z 2 v = z^z^Zq 1 , 

V-'ZsV = Z^Z^ 1 , = Z , V^V = Z^ 1 , V" 1 ^ = Z 6 - x , (21) 

which can be further simplified by using fll8|) . In the additive form ( [H|)and (|19"D can we 
written in full generality as: 

6 

h- 1 z t h = Y,Mh)z k , fceff (22) 

fc=0 

where .P(/i) is a representation mod 2 of H. 
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Moreover the group H has an external automorphism v given by : 

u(y) = v- 1 = -v 

u(U) = U- 1 

v{Z p ) = Zll v = Z^ v (23) 

V. ORBITS OF HOMOLOGY AND COHOMOLOGY GROUPS 

As anticipated, the analysis of orbits and of external automorphisms plays a central role 
in the computation of Z d . We thus give here the explicit construction of such orbits in Hi 
and together with their duality map. 

A functional G H 1 can be defined by the equivalent conditions: 

6 

a(gZi) = eia(g) ,(f)(Zi) = ei , i = 0,...,6 ; JJ e< = 1 (24) 

1=0 

is then identified by (e , . . . , e 6 ) . Let ( G Hi and let h^Ch be the translated cycle. Then 
the translated functional <fih is defined by 

( f )h (h- 1 Ch)=<j ) (0 . (25) 



We have then: 



MU-'Z.U) = <j>{Zi) =e l = MZi+i) , (26) 



and hence 



4>u(Zi) = ej_i , fiu^Zi) = e i+ i . (27) 



Similarly we find 



0y(\/- 1 Zo\/)=0(Zo)=eo = 0y(Z 4 ) , (28) 



and 
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(f)v{Zo) = e 4 , (fiv(Z 3 ) = e x e 4 , (fiv(Z 2 ) = e e 4 e 2 , 



4>v(Zi) = e 3 e , (pv(Z 5 ) = e 5 , (fiv(Z 6 ) = e 6 



(29) 



There are 2 2g = 64 (g = 3 denoting the genus of A) different elements in H 1 which fall 
in 5 different orbits containing 1,7,7,21,28 elements . We list these orbits with only one 
signature out of each cyclically permuted septet: 



A : (1, 1, 1, 1, 1, 1, 1) (trivial) 
B: (1,1,-1,1,-1,-1,-1) 
C: (-1,-1,-1,1,-1,1,1) 

D : (-1, -1, -1, 1, -1, -1, -1) , (-1, -1, 1, 1, 1, -1, -1) , (1, 1, -1, 1, -1, 1, 1) 
E : (-1, 1, -1, 1, -1, 1, -1), (1, -1, -1, 1, -1, -1, 1), (-1, 1, 1, 1, 1, 1, -1), 

(1,-1,1,1,1,-1,1) (30) 



The orbits B, C are mapped into each other by the action of v. This proves incidentally 
that v is external. We also list in similar fashion the dual orbits in Hi : 



A : (0, 0, 0, 0, 0, 0, 0) (trivial) 
B: (1,1,0,0,1,0,1) 
C: (1,0,1,0,0,1,1) 

D : (0, 0, 1, 0, 1, 0, 0) , (0, 1, 1, 0, 1, 1, 0) , (0, 1, 0, 0, 0, 1, 0) 
E: (1,0,0,0,0,0,1), (1,1,1,0,1,1,1), (1,1,0,0,0,1,1), 

(1,0,1,0,1,0,1) (31) 



as it can be checked by using (p8|)-(p9|). 



For generic elements £, (' G H\ we use the additive form: 



(32) 



i=0 



i=0 



and define the intersection number (mod2) 
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r(C,C')=T(C',C) = (-l) n(C,C0 • (33) 
Q is the Z 2 valued form of eq.(|5]) explicitly given by: 

«(C,0 = E x(*-*0"i™* (34) 

i,fe=0 

with x(i) = 1 if % = 1, 2, 5, 6 (mod7) and x(i) = otherwise, as it can be verified on the 



graph A*. r((,(') is invariant under conjugations by H, see fl2lD , and is multiplicative, i.e. 

t((,C + (")=t((,C)t((,(") • (35) 

We define also: 

0(() = (_l)3Ei,*=oX(i-fc)nin* ^ (36) 
so that #(Zj) = 0(1) = 1 and 9(Q is a sign functional as 77 in (j5|): 



0^ + = 0(00(0 (CO • (37) 

All these definitions can be extended naturally to a wide class of lattice groups. 
Let then : 

6 

m i = E X(i ~ k ) n k (38) 

fc=0 

and define the duality map ip ,Hi —> H 1 : 

^(l>^) =((-l) mo ,.-, (-l) ms ) (39) 

The 2 29 elements of H^H 1 are labeled as ,where (p(Ci) — <Pi , by an index / = 
1, ...,2 2g with (i,4>i the trivial elements and sorted in such a way that in (|3"T|) , (|3H|) (p maps 
corresponding orbits in Hi,H x . Clearly then 0/ (( K ) = t(Q, (k)- We set T]i (() = <pi (() 9 (£) 
, si = s Vl so that we may label Pfaffians equally well with elements G H l and rewrite (pTj|) 
as: 

Z d = 2- g '£s I Pf((f> I ) . (40) 
14 



The importance of orbits should be clear once we realize that in fl40D for elements (pi, <fij 
in the same orbit we have Pf(<pj) = Pf(<f)j). In this way the effective number of different 
terms in Z& reduces to the number of orbits in H 1 . 

Further reductions arise from external automorphisms. When omitted, as in (|10D , we 
assume summation ranges on /, K to be 1...2 29 . Since ip commutes with the group operations 
it maps orbits in H\ into dual orbits in H 1 and we use the same label A...E for pairs of dual 
orbits. Moreover the 2 2a x 2 2a matrix $ = $ T of elements <&ki = 2~ 9 t((j, (k) is orthogonal. 
We have first of all: 

E^/ = 2- 2 *£i = i . (41) 

We have next: 

p(k, k') = E $ki*k>i = 2- 29 E MCi)MO) = 2" 29 E MO + 0<M0 + (42) 
11 1 

1 

therefore p(K,K') = if 3f : <M0<M0 ± 1- But if : <M0<M0 = 1 th en K = K' 
therefore p(K, K') = if K ^ K' . This proves that $ is orthogonal and that $ = 

We conclude this section with some rather technical formulas whose role will be crucial 
in what follows (sect. VIII). 

The sum J2k 0(Ck) can be readily evaluated by using a standard basis for Hi given by 
Xi, Yi , i = l...g such that 

0(E (PiXi + q l Y i ) ) = (-l)^tx (43) 

i 

where for instance: 

X l = Z Q ,Y X = Z U X 2 = Z x + Z 5 ,Y 2 = Z + Z 2 + Z 5 ,X 3 = Z + Z 3 ,Y 3 = Z 1 + Z 3 + Z 6 (44) 

In this form J2k 9(Ck) factors into g independent and equal sums each yielding a factor 2 
and hence Ek&(Ck) = 2 9 ■ 
Furthermore: 
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2" 9 E <Ci, CkMCk) = I' 9 6(0) E 6 ((k + 0) = 0(0) ■ (45) 

K K 

since (k + Ci runs over the whole group H\ taking every element once just as (k and therefore 

E*0(Cjr + Ci)=E*:0(C*) = 2 fl . 

VI. IRREPS OF H 

The group H has 168 x 64 = 10752 elements and in order to perform harmonic analysis 
and obtain partial block diagonalization of X(4>) and of Pfaffians we must find the unitary 
irreducible representations (irreps) of H. 

The trivial functional (1,1,1,1,1,1,1) must be dealt with separately and requires the 
construction of the irreps of the factor group H of H and central extension of G defined by 

U 7 = Q, V 2 = Q, (UV) 3 = Q, (y[/ 3 ) 4 = Q, (46) 

(i.e. Zi = 1) of order 168 x 2 = 336. Therefore we have the relators : 

a(gZi) = a(g) , a(gQ) = -a(g) . (47) 

Because of (f4~7|), only a subset of the unitary irreps of Hq is actually used in the harmonic 
analysis. In order to see it let us write such irreps as D^(g) with J a convenient label and 
a,/3 = 1, . . .,dj where dj is the dimension of the irrep . The matrix elements satisfy the 
orthogonality relations: 

E DipigfD^^g) = 5 JJ '5 aa ,5 m . (48) 

Let us define 

aj> = E <9)D J aP (g) . (49) 

Applying fl49|) to fl47|) we get 

E <gQ)D J ap {g) = -< = a^D^Q) . (50) 
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But Q is a central element and Q 2 = 1. Therefore by Schur lemma D^(Q) must be 
proportional to the identity, i.e. 

D^(Q) = k5 j(3 ,k = ±1 , (51) 

therefore from fl50f ) 

^< = -< , (52) 

i.e. k = — 1. It follows that only irreps having k = — 1 actually contribute to the Fourier 
expansion of a(g). Those irreps which are also irreps of the original group G (i.e. those with 
K — 1) are absent from the expansion. There are only 5 irreps of H with k — — 1 of dimension 
4, 4, 6, 6, 8, satisfying separately the Burnside condition 4 2 + 4 2 + 6 2 + 6 2 + 8 2 = 168. The 
detailed form is listed in Appendix A. The matrix elements of these irreps are polynomials 
in K = exp(iir/7) with K 7 — — 1. By abuse of language we may write 

Q = -l, U 7 = -1 , V 2 = -1 , {UVf = -1 , (VU 3 ^ = -1 , (53) 

instead of (|46|). 

As for the remaining irreps of H the most convenient way is to obtain them as induced 
representations on the cosets of the subgroup L of H generated by: 

V , U- l VU 2 , U- 6 VU 3 , U-*VU A , U- 5 VU 5 , (54) 

as well as the cosets of the group v{L) generated by : 

V~\ UV^U- 2 , U^V^U' 3 , U 4 V- l U-\ U^V^U- 5 . (55) 

All generators and their inverses are of the form U a ^VU b , b = 0, . . . , 6. Let set briefly 
v — V and t = U^VU 2 ; we find then 

U-*VU 3 = vt , U~*VU A = Z 1 (tv) 2 = (tv) 2 Z 4 
U- 5 VU 5 =vt(vt~ 1 ) 2 Z 1 = Z 2 Z 6 vt(v r 1 ) 2 

(vt) 4 = QZ 3 , {tv) A = QZ X Z± . (56) 
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It can be verified that {v, t, Zi, Q} generate the whole subgroup L. Let suppose that a 
representation A : L —>■ Hom(£) is given on a linear space L of dimension n. Then A can be 
extended to the induced representation /x : H — > Hom(7i) where 7i = Q) P=Q £ P and C p are 
isomorphic to C Let ipi be a basis on C, ijj^ p a basis on C p and ipi, p +7 = —ipi, P ■ The action 
of [/ on Ti is then defined by 

Uip itP = ipi, p+ \ , (57) 

hence U 7 = — 1. Moreover we set 

W^,o = ^k,ph,ii w ) > w E L . (58) 
From (|5B|)-(|5BD derive the representation of V on 7i as follows: 



^V fc ,oA(f/- a(6 Vf/ 6 ) fc , l = ^ M 6)A(£r°<>W) fc)< . (59) 

Since b can take all values 0, . . . , 6 we deduce the complete representation of V on TC. The 
problem reduces now to that of finding a suitable representation of L on C. Should we set 
Q — Zi — 1, L becomes the octahedral group T(2, 3,4) of order 24. We are interested in 
irreps of L in which all \{Zi) are diagonal and X(Zf) = 1. Each X(Zi) is then a diagonal 
block matrix. In fact we have 

Ztflnj = Z U p ^ lfi = U p Z p ^ lfi = U p X(Z p )^ fi = X(Z p )^ ijP . (60) 

It is then clear that /x( Z s ) is the diagonal matrix with entries X(Z S ), X(Z s+ x),. . ., A(Z s+6 ) 
, where X(Z S ) has n eigenvalues = ±1, i — 1, . . . , n , s = 0, . . . , 6. It follows also that 
each vector ip itP has eigenvalues 

Z s ipi :P = £,i jS +pipi,p > (61) 
which define the element 5 e H l : 

^ = (£i,p> • • • j £i,p+6) (62) 



For the group v{L) the corresponding formulae are obtained by setting v' = V 1 , t' — 
UV^U- 2 . We have then 

U 6 V~ l U- 3 = v't' , U A V~ l U~ A = Zz {t'v'f = {t'v'f Z Q , (63) 
{v't'f = Q Z~\ {t'v'f = Q Z, l Z, 1 . (64) 

If an irrep of H contains a signature H it contains also the whole orbit including its 
cyclical permutations. Non trivial functionals have 7 distinct cyclical permutation of any 
H corresponding to p — 0, . . . , 6 in fl62|) . The irrep fx oi H generated by A has dimension 
In and as many signatures. If it contains different orbits then is reducible. Therefore 
irreps can be grouped in disjoint subsets characterized by orbits. Irreps belonging to orbits 
A, D, E can be obtained either from Lor v{L) whereas B and C can be obtained only from 
L and v(L) respectively. We list in appendix B the explicit representations A of the group 
L corresponding to each orbit by giving \{y) , A(t) and \{Z{). In this way we obtain the 
irreps of the orbits A, B, D, E. Irreps of orbit C can be obtained by applying v to irreps 
of B. Clearly orbits may appear in some irreps with multiplicities m = 1,2,4 where n/m 
is the number of inequivalent signatures under cyclical permutations contained in the orbit. 
The method applied to orbits A yields representations which reduce into the irreps already 
discussed. The dimensions of the irreps satisfy a separate Burnside condition 

E = <0)N (65) 
Jeo 

where n(0) is the number of signatures belonging to a given orbit O. 

VII. THE DECORATED LATTICE AND THE PARTITION FUNCTION 

We need now a more explicit description of the decorated lattice A* along lines already 
discussed in ref. JJ. Each element h G H identifies uniquely a site of A which we also label 
with /i,being intended that hQ ~ hZi ~ Zih identify the same site of A. h is connected 
to other 3 sites V h,Uh,U~ 1 h and is replaced in A* by 3 sites h a ,hj),h c to which we as- 
sociate Grasmann variables a(h) , b(h) , c(h) . The original oriented bonds h — > Vh,h — > Uh 
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are replaced by oriented bonds h a — > Vh a ,hb — > C//i c (with exchange interaction Jy,Ju 
respectively) together with the additional decorating bonds h a — > h c , h c — > /i 6 , ft, 6 — > /i a . For 
each Pfaffian we need boundary conditions of the kind discussed above: 

a(hZi) = 6ia(h) , a(Qh) = a(hQ) = —a(h) (66) 

and the corresponding conditions obtained by replacing a with b, c. Under these assumptions 
the 2-form @ can be rewritten explicitly as : 

f =J2 ( a ( h ) A c ( h ) + c ( h ) A K h ) + K h ) A a(h) + | a(h) A a(Vh) + x b(h) A c(Uh)) (67) 

where y = coth(/3Jy), x = coth(/3J[/) . 

Harmonic analysis on A can be performed by recalling the Fourier components of 
a(h) , b(h) , c(h) given by a(h) = Y^aB^aB{^)* a aBi e ^ c - • To this purpose we consider the 
matrices a J ,b J ,c J of elements , b^ g , c£g where a, (3 = l...dj . By using the orthonor- 
mality relations ( Pf ) satisfied by the matrices D^{h) we rewrite / = J2j f J where 

f J = Tr (a Jt A c J + c Jt A b J + b Jt A a J + | a Jt A D J (V) a J + x b Jt A D J (U) c J ) (68) 

and Jt labels the complex conjugate irrep of J. The key point is that different pairs J, Jt 
lead to disjoint sets of Grassman variables, moreover f J + f J ^ separates into the sum of 
partial forms Y^,a{f J )aa + c.c. . Therefore the final Pfaffian is the product of dj partial (and 
identical) Pfaffians, which must be computed explicitly by use of (|8]). Once this is done we 
must determine the coefficients si appearing in ([40D. Since all the sj and Pfaffians Pf(4>i) 
sharing the same orbit O of (pi are equal we can set So = n(0) sj , Pf(4>i) = pf(0) , V0j G O 
and P0| ) reduces to: 

Z d = 2-^SoPf(0) . (69) 
o 

The sign in Pf(<fii) of a term associated with the path ( is given by 777(C) = t(Cj> C) ^(0- 
The si must be chosen in such a way as to set the final coefficient of C,k in the expansion 
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2- s $>t(C/,Ck)0(Ck) (70) 
equal to 1. Taking into account ([45]) we see that the condition 

V/ : 2~" E si r(Ci, Ck) = E s i = 0((k) (71) 
i i 

has the simple solution: 

= e(0) , (72) 

and thus S A =1,S B = S = 7, S D = 21, S B = -28 and S A + S B + S c + S D + S E = 8 = 2». 

VIII. PFAFFIANS AND HARMONIC ANALYSIS ON THE LATTICE 

The boundary conditions (|66|) on the Grassmann variables a(h) , b(h) , c(h) imply strong 
restrictions on the Fourier components a J af3 , b J a[3 , c J a[3 similar to those already derived for 
Q . Consider now Pf (0) where (ft = (eo, ei, £2, €3, €4, £5, eg). From ( |66D we get : 

E a(hZ t )D J ^(h) = e t a J ai = a^D^Z,) (73) 

In all the explicit representation of H listed in the Appendix B the matrices Di^iZi) are 
diagonal and 

Dt P {Z i ) = 4(>y)5tf . (74) 

Clearly a^ 7 ,6^ 7 ,c^ 7 = unless ef(j) = ej , z = 0...6. The orbit O of identifies a subset 
of irreps and in (|68"D only a rectangular submatrix of a J , b J , c J of dimension c?j x dj /n(0) 
survives where dj is the dimension of the irrep and dj /n(0) the number of times a given 
signature in O appears in the irrep. The total number of surviving components of a J a ^ 
among all irreps is in any case equal to the order N of the group G, i.e. 168. 

Each a J ,b J ,c J is therefore partitioned into dj /n(0) rectangular submatrices labeled 
by the signatures in O. Only one of these matrices needs to be used because they all yield 
the same final expression for the Pfaffian. Since in (|6~8"D all partial forms fj a yield identical 
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Pfafnans we may consider just one contribution and drop the index a. In this case the 
square of the partial Pfaffian can be written as the determinant of the block matrix Aj 



A. 



yD J {V) 1 -1 

-1 l-xD J {U) 

1 -l+xD J (U~ l ) 



(75) 



Notice that Aj = —Aj. Upon multiplication of the last group of columns by ( 1— x D (U) 



and the last group of rows by (l— xD J {U 1 



-i 



we find as in 



Det (Aj) = Det [y (l- x D J (U)) D J (V) (l- x D J '(U' 1 )) + x (D J (U) - D J (U~ 1 ))] (76) 

and therefore the effective dimension of the final determinant reduces to dj. The dimer 
generating function Z& is then: 



Z d = 2S J2 So Pf(0) = 2-° J2So II p f. 



dj/n(0) 



(77) 



The degree of Pfj in x, y is dj,dj/2 and therefore the total degree of pf(0) in x is 
J2jeo dj/n(0) = N because of the Burnside condition and half as much for y. Not all 
Det(Aj) are distinct, besides Det(Aj) = Det(Ajt) the external automorphism v maps 
pf{B) into pf(C) and viceversa thus pf(C) = pf(B) . The number of independent Pfafnans 
is therefore further reduced to four (64 — > 5 — > 4) by the existence of external automor- 
phisms. 



If dj j n(0) is odd the sign of Pfj = ±J Det (Aj) must be fixed unambiguously by 
computing it for instance in the limit f3 — *■ 0. becomes then a polynomial in x, y of 
degree N,N/2 [fT2| . In Fig.2 we give the plot of the specific heat corresponding to a ferro- 
magnetic choice of the exchange interactions, (Ju = Jy = 1). The specific heat presents a 
very sharp peak which prefigures the transition form and ordered ferromagnetic phase to a 
disordered paramagnetic one. When Ju, Jy have different signs or are both antiferro mag- 
netic (i.e. negative) the model is totally frustrated and the specific heat becomes smoother. 
Moreover, different finite values (> log2/iV) of the ground state entropy are found, signaling 
an exponential degeneracy of the ground state. 
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IX. CONCLUSIONS 



In this paper we have developed a formalism capable of dealing efficiently with Ising 
models defined on group lattices of non trivial genus g. In particular we have applied the 
method to the Ising model on the Klein group L(2, 7) having g = 3 and N = 168. We found 
that the computation of the partition function Z is greatly simplified by use of symmetries 
of an extended group, both internal and external to the group, which reduce the number 
of and provide explicit formulas and topological interpretation for the sign and weight of 
Pfaffians in the expansion of Z. We plan to apply this method to other lattices where N is 
large and g is comparable to N. 
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Fermionic irreps of Hq : 



= diag(-l,K,-K 2 ,-^ 4 ), = ([/«)* , 
= diag(K, -K\ -K\ -K\ K\ K 3 ) , = (u^)' 
U {5) = diag(tf, -K 2 , -K\ -1, -is: 6 , if 5 , K 3 , -1) . 



V 0) 



2i 



1/2 -i- -i- 



Cl c 3 c 2 



75 C3 C2 Cl 
V 73 C2 Cl C3 



with c fc = cos (^p) ,Jfe = 1,2,3. V< 2 > = (V« 



V® = 



2i 



A B 

B -A 



where 



^ S 2 - Si - S 3 - V^Sl ^ 



B = r 



Si - S3 - V^Sl S 2 - V2-Ss 

s 3 - y/2s 1 s 2 - V%s 3 Si - a/2s 2 
/ 

s 3 + (y/2 - l)si s 2 + (v^2 - l)s 3 si + V2 - l)s 2 

s 2 + ( v / 2-l)s 3 s 1 + v / 2-l)s 2 s 3 + (v^-l)si 

y si + - l)s 2 s 3 + (v^ - s 2 + (v^2 - l)s 3 



\ 



with si = sin (f ) , s 2 = - sin (^f) , s 3 = - sin (^) and r = y 2 + ^2 . 
Finally V< 4 > = (V( 3 ))*and 



\/( 5 ) = - 
7 



C D 

D -C 



where 
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c = 



si + s 2 

51 + s 3 

5 2 + S3 



|(s 2 - Si) 



51 + S 3 

5 2 + s 3 
Sl + s 2 



|(si - s 3 ) 



s 2 + s 3 
Sl + s 2 



|(s 2 -si) ^ 
3/ si - s 3 ) 



2 1 

si + s 3 y / |(s 3 -s 2 ) 

(s 3 - si) |(si - 2s 2 + s 3 ) 



and 



s 2 


Sl 


S3 




Sl 


S 3 


s 2 


7l( S 3-Si) 


S 3 


s 2 


Sl 


75( s 2 - s 3 ) 



D = 2V3 

v 7 2 = (si-s 2 ) 73(53 -s x ) ^(s 2 -s 3 ) |(si - s 3 ) j 

All these irreps are changed into equivalent and conjugate expressions by replacing K 
with —K 2 . 

Appendix B 

Putting \{x) = xo for the representations A belonging to the orbit O we have 
Case B: 

( (l+i V2) _i \ 





( . 

% 




1 


V B = 














\ 



(l^jvg) 
2 



Case C: 
Case D: 



V C = Vg 



t c = t 



-1 

B 1 



Z\r — Z 



l,B 



V D 



1 

72 



'00-11 



0-1-10 

1 1 
-1 1 

-% -% 
-% % 



£ (1) 



1 

72 









V2 




















V2 




















— i 


— i 














—i 


i 


—i 


—i 














—i 


i 















Z«=diag(l, 1,1, 1,-1,-1), 
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,(2) _ / ,(2) _ / 7 (2) 7 (1) . 



Case : 



7 (i) 



/ -i N 
-i 
-i 
^ -i 
diag(l, -1,1,-1) 



t (i) 

"E 



,( 2 ) 



^-1 0^ 



-% 
1 
-i 



(4'T 



(2) 



(4 1J )* 



7 (2) 



7 (1) 



.(3) 

1,B - 



f 











1 





0^ 




f 


1 

2 


-V3 
2 

















\ 

















1 









V3 
2 


1 

2 


























i 






























1 








o 











—i 











/(3) _ 
, l E — 



















1 






































-V3 
2 


-1 










































-1 

2 





-1 







































2 


1 

2 

















i 















1 

2 


2 














V 





















\ 








2 


-1 

2 











1 


diag(- 


-1,- 


-1,1,1,- 


-1,- 


1,1,1) . 
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FIG. 1. Intersection between I and £' and schematic representation of (,(' an d C + C ( or 
and £"). 



FIG. 2. Specific heat versus temperature for J\j = Jy = 1. 
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